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Abstract. The understanding of present atmospheric trans- 
port processes from Southern Hemisphere (SH) landmasses 
to Antarctica can improve the interpretation of stratigraphic 
data in Antarctic ice cores. In addition, long range trans- 
port can deliver key nutrients normally not available to ma- 
rine ecosystems in the Southern Ocean and may trigger or 
enhance primary productivity. However, there is a dearth of 
observational based studies of dust transport in the SH. 

This work aims to improve current understanding of dust 
transport in the SH by showing a characterization of two 
dust events originating in the Patagonia desert (south end of 
South America). The approach is based on a combined and 
complementary use of satellite retrievals (detectors MISR, 
MODIS, GLAS , POLDER, OMI,), transport model simula- 
tion (HYSPLIT) and surface observations near the sources 
and aerosol measurements in Antarctica (Neumayer and 
Concordia sites). 

Satellite imagery and visibility observations confirm dust 
emission in a stretch of dry lakes along the coast of the 
Tierra del Fuego (TdF) island (~54° S) and from the shores 
of the Colihue Huapi lake in Central Patagonia (~46° S) in 
February 2005. Model simulations initialized by these ob- 
servations reproduce the timing of an observed increase in 
dust concentration at the Concordia Station and some of the 
observed increases in atmospheric aerosol absorption (here 
used as a dust proxy) in the Neumayer station. The TdF 
sources were the largest contributors of dust at both sites. 
The transit times from TdF to the Neumayer and Concor- 
dia sites are 6-7 and 9-10 days respectively. Lidar obser- 
vations and model outputs coincide in placing most of the 
dust cloud in the boundary layer and suggest significant de- 
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position over the ocean immediately downwind. Boundary 
layer dust was detected as far as 1800 km from the source 
and ~800 km north of the South Georgia Island over the cen- 
tral sub- Antarctic Atlantic Ocean. Although the analysis sug- 
gests the presence of dust at ~ 1500 km SW of South Africa 
five days after, the limited capabilities of existing satellite 
platforms to differentiate between aerosol types do not per- 
mit a definitive conclusion. In addition, the model simula- 
tions show dust lifting to the free troposphere as it travels 
south but it could not be confirmed by the satellite observa- 
tions due to cloudiness. 

This work demonstrates that complementary information 
from existing transport models, satellite and surface data can 
yield a consistent picture of the dust transport from the Patag- 
onia desert to Antarctica. It also illustrates the limitation of 
using any of these approaches individually to characterize the 
transport of dust in a heavily cloudy area. 


1 Introduction 

Airborne dust can impact significantly the radiative balance 
of the Earth’s atmosphere (Tegen et al., 1996; Sokolik et al., 
2001; Miller et al., 2004). In addition it can impact marine 
and land ecosystems. For example, a number of studies have 

shown the effect of long range transport of nutrients carried 
by dust to the Amazon forest (Swap et al., 1992; Okin et 
al, 2004) as well the disruption of primary productivity in 
the North Atlantic (Mallet et al., 2009) by reducing solar 
radiation reaching the ocean surface. Long range transport 
of dust can supply key nutrients (Jickells et al., 2005) that 

stimulate primary productivity in areas where phytoplankton 
growth may otherwise be minimal. This is relevant in the cur- 
rent climate change context because oceans are responding 
to global warming by becoming more stratified in the upper 
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layers. This may prevent the surface ocean from mixing with 
the deep ocean, thereby decreasing the nutrient fluxes from 
below (Levitus et al., 2000; Sarmiento and Gruber, 2006) and 
affecting oceanic primary production and efficiency of CO 2 
transfer (Behrenfeld et al., 2006). This climate- warming in- 
duced situation makes oceanic primary production more de- 
pendent on the nutrient input from external sources, such as 
atmospheric deposition. This fact makes the study of the 
presence of dust over the Southern Ocean (40S to 60S) of 
particular interest because its role in the carbon cycle (it con- 
tains some of the largest CO 2 sinks, Takahshi et al., 2009) 
and the concerns on how it would respond to a change in 
nutrient availability (Watson et al., 2000). 

The hypothesis of dust inducing an increase in primary 
productivity was proposed by Martin (1990) as an explana- 
tion of the observed decrease of atmospheric CO 2 with an 
increase deposition of dust found in ice cores at the end of 
the cold periods (Petit et al., 1999). In addition, dust layers 
in ice cores are one of the few proxies available to study paleo 
atmospheric transport phenomena (Fischer et al., 2007; Lam- 
bert et al., 2008). Isotopic analysis of ice cores show that the 
Patagonia desert was amongst the largest contributors to the 
dust found in East Antarctica (Basile et al., 1997; Delmonte 
et al., 2009). Although dust deposition was more abundant in 
glacial times, recent studies in snow samples in the Antarctic 
Peninsula (McConnell et al., 2007) and East Antarctica (Sala 
et al., 2008; Marino et al., 2008a; Lambert et al., 2008) indi- 
cate there has been recent deposition of dust (in the last 100 
years) and it was suggested that the Patagonia desert was an 
important source along with Australia (Revel-Rolland et al., 
2006; Marino et al., 2008a). 

Much of what is known about transport of tracers in the 
high latitude Southern Hemisphere is known through model 
studies (Jickells et al., 2005; Li et al., 2008; Mahowald et 
al., 2009). Observational studies of aerosol transport into the 
Southern Ocean have been notably missing and restricted to 
occasional research cruises (Baker et al., 2006; Wegener et 
al., 2008). The high latitude marine environment is rather 
cloudy and dust emissions from land masses in the SH are 
sporadic and not very abundant. Both features contribute to 
a low probability of detection with satellite detectors because 
detection algorithms require clear sky pixels for dust detec- 
tion. In addition, there are very few data sets of continuous 
aerosol monitoring in the Southern ocean region and Antarc- 
tica. Only recently and thanks to the increase of in- situ mon- 
itoring (Gaiero et al., 2003; Gaiero, 2007) and the overabun- 
dance of satellites with complementing detection capabili- 
ties, a number of combined observational and modeling stud- 
ies have been published (Gasso and Stein, 2007; Fiebig et al., 
2009). Also, detailed studies of modern aerosol composition 
in coastal and central Antarctic regions have become avail- 
able recently (Udisti et al., 2004; Fattori et al., 2005; Becagli 
et al., 2009) 

This paper presents a detailed study of two dust events 
originating in the south end of Patagonia. Using a combina- 


tion of model and satellite data analysis, this work tracks and 
characterizes the progression of a dust cloud over the South- 
ern Ocean and compares the model outputs with aerosol mea- 
surements at two sites in East Antarctica. 

This work is motivated by the direct observation of a sharp 
peak of non-sea-salt Ca 2+ (nssCa, a proxy for dust, Fischer 
et al., 2007) deposition during a single day (9 March 2005) 
at the Concordia Station (75.1 S, 123. 35E; 3233ma.s.l., 
Antarctica), the location of the EDC deep ice core, drilled in 
the framework of the EPIC A (European Project for Ice Cor- 
ing in Antarctica). The short period of exposure permitted 
the exploration of the provenance of the dust using a back- 
trajectory modeling approach and a detailed look at satellite 
imagery. 

This paper is structured as follows. Section 2 details the 
modeling tools and data sets available and their main features 
and limitations, section 3 explains the approach used to an- 
alyze the data as well as a description of model outputs and 
observations. Section 4 summarizes and integrates the find- 
ing of Sect. 3. Section 5 discusses the main conclusions and 
implications of these findings. 

2 Data 

2.1 Satellite data 

The satellite data from a number of detectors was selected for 
aerosol type differentiation, dust event localization and 3-D 
tracking after the event. The satellite data analyzed was from 
the detectors MODIS (Remer et al., 2005), GLAS (Spinhirne 
et al., 2005), MISR (Diner et al., 2001), PARASOL (Herman 
et al., 2005) and OMI (Torres et al., 2007). The main capa- 
bilities of these detectors have been documented elsewhere 
and only those features relevant for this analysis will be dis- 
cussed. All satellite data used in this study were extracted 
from the standard products (visible images and level 2) made 
available for each satellite. 

One important hypothesis in this study concerns the detec- 
tion of aerosol type present. The high latitude southern hemi- 
sphere marine environment is mostly dominated by marine 
aerosols (largely derived DMS and sea salt), with occasional 
seasonal intrusions of dust and smoke from biomass burn- 
ing (Kaufman et al., 2002). Although pollution is important 
in some major metropolitan areas in the mid-latitudes, these 
aerosols are not known for traveling long distances in the 
Southern Hemisphere. Large scale transport of smoke usu- 
ally occurs in late August in West Africa and August to Oc- 
tober in the Amazons with the latter is known to reach high 
latitudes (>60 S) occasionally (Pereira et al., 2006; Fiebig et 
al., 2009). During the period of interest (end of austral sum- 
mer), dust is a sporadic feature in the southern hemisphere 
deserts and biomass burning is uncommon. 

From the viewpoint of remote sensing, there are specific 
satellite markers that can be used to differentiate dust from 
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sea-salt. For example, the Aerosol Index (AI) product from 
OMI (Torres et al., 2007) is sensitive to absorbing aerosols, a 
common characteristic of dust and not of sea-salt aerosol. 
MODIS fine mode fraction aerosol product (a measure of 
which aerosol mode is dominant) can be used as a marker 
for dust provided that high aerosol concentrations are present 
(Kaufman et al., 2002). The MISR non-spherical fraction 
(Kalishnikova and Khan, 2008) provides information on the 
fraction of aerosol optical depth that can be matched to a non- 
spherical aerosol model. Likewise, the POLDER aerosol al- 
gorithm produces a dust proxy called non-spherical fraction 
too (Herman et al., 2005). 

However, all these proxies have limitations that become 
apparent in the southern environments. Abundant cloudi- 
ness and the resulting contamination is a major limitation 
that results in very few clear sky pixels. In this case, OMI 
and PARASOL retrievals are particularly affected. Gener- 
ally speaking, retrievals near the source (usually on the day of 
emission and the day after) are robust if there is enough cloud 
clearing and aerosol concentrations are high. It is expected 
that the five sensors would be consistent in observing the 
same aerosol although they may not detect it as dust. How- 
ever, the tracking of the dust as it advects is a challenging 
task. Dispersion and dilution of the dust cloud results in low 
aerosol concentrations and very low reflectances reaching the 
satellite sensor. Further some of the aerosol type proxies re- 
trieved have not been compared or established their range of 
validity by comparing with independent measurements (AI, 
non-spherical parameters from PARASOL and MISR). Other 
like the MODIS fine mode fraction and Aerosol Index prod- 
uct has been widely used and compared against observations 
(Kleidman et al., 2005). 

All satellite data products were used at level 2, that is, 
non-gridded native parameter product except in the case of 
PARASOL. The data analyzed was processed by the respec- 
tive science teams using the following data algorithm ver- 
sions: MODIS: collection 5, MISR: version 12 (F12_0022), 
OMI: collection 3, PARASOL: Gridded atmospheric prod- 
uct, software version 10.5, GLAS: version 2.8. 

2.2 Surface data 

Different sources of surface data were used in this analysis. 
Model outputs were compared with surface measurements 
of aerosols at two sites in Antarctica. In addition, the model 
simulations were initialized and constrained by meteorologi- 
cal information from several weather stations along the coast 
of Patagonia. These data sets are described next. 

2.2.1 Concordia Station in Antarctica 

The Concordia Station (75.1 S, 123. 35E; 3233ma.s.l., 
Antarctica) is part of the EPIC A (European Project for Ice 
Coring in Antarctica) project and it is where Dome C deep 
ice core has been drilled. A long term aerosol monitoring 


program is active since 2004. The program consists on the 
collection of aerosol samples at different temporal resolu- 
tion, with and without particle size separation, in order to 
obtain highly detailed information on aerosol chemical and 
physical properties, as well as changes in transport and de- 
position patterns in different seasons (Becagli et al., 2009). 
The study of mechanisms of present-day aerosols transport 
and deposition is basic for the correct interpretation of ice- 
core proxies, among which dust constitutes one of the most 
studied (Lambert et al., 2008). Rapid events of dust trans- 
port to Dome C are studied by measuring nssCa as soluble 
proxy of dust in aerosol samples collected at higher tempo- 
ral resolution (from daily to weekly). NssCa 2+ concentra- 
tions in ice cores are commonly used to infer information 
about terrestrial inputs in polar regions (Wolff et al., 2006), 
even if the quantitative relationship between nssCa 2+ and 
dust are not fully defined for different climatic regimes, de- 
pending on dust mineralogy (e.g. Bigler et al., 2006; Ruth 
et al., 2008); however nssCa 2+ is calculated from Na + and 
Ca 2+ concentrations (using the formula nssCa 2+ = Ca 2+ - 
(Ca 2+ /Na + )sw * Na + , where sw stays for sea water), mea- 
sured with high sensitivity in ice and aerosol samples by Ion 
Chromatography (Morganti et al., 2007). Data presented in 
this work come from nssCa 2+ measurements performed on 
aerosol samples collected by 8-stage Andersen impactor for 
size-resolved (from 0.4 to 10 pm) aerosol collection. Filters 
are exposed every one or two weeks using a pumping volume 
of about 28.3 1 min -1 . Data plotted in Fig. 1 refer to the sum 
of nssCa 2+ concentrations in the 8 stages, i.e. they constitute 
the atmospheric concentration of nssCa 2+ soluble fraction of 
particulate in the range 0.4-10 pm. In the particular case of 
the main event of dust deposition studied here, aerosol col- 
lection started at about 3 p.m. (Dome Concordia local time, 
i.e. UTC +7) of 8 March, and ended approximately the day 
after (9 March) at the same time. On 9 March 2005, operators 
had to stop sampling operations temporarily for instrument 
maintenance and removed filters which had been installed 
just 24 h before. The sampling resumed normally after a few 
days. During the post campaign analysis, an anomalous peak 
in nssCa was detected in all the filters exposed in that 24hrs 
period (largest concentration seen in Fig. 1). 

2.2.2 Neumayer surface station in Antarctica 

The Alfred Wegener Institute for Polar and Marine Research 
(AWI) runs a permanently manned Georg-von-Neumayer- 
Station (70.65 S, 8.25 W, 36ma.s.l) at the northwest edge of 
the Weddell Sea in Antarctica. It is part of the WMO Global 
Watch Stations (Konig-Langlo and Weller, 2006). Aerosol 
sampling was made at the Air Chemistry Observatory, about 
1.5 km south of Neumayer station. Weller et al. (2008a) 
reported surface measurements of aerosol number concen- 
tration, optical properties (scattering and absorption coeffi- 
cients, angstrom exponent) and limited chemical composi- 
tion. The data was reported as daily averages from January, 
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Fig. 1. Non sea-salt calcium extracted from airborne aerosols mea- 
sured at the Concordia Station site (Dome C, Concordia project) 
during the 2004-2005 campaign. The largest peak corresponds 
to a 24 h filter measurement from 09:00 UTC 8 March 2005 to 
08:00 UTC 9 March 2005. 


2004 to December, 2006. The version of data used in this 
study is the one available at the free public archive Pangea 
(Weller et al., 2008b). Because no specific chemical anal- 
yses for dust were available at Neumayer, the aerosol ab- 
sorption coefficient can be used as proxy of dust if making 
reasonable assumptions. Dust is more absorbing than marine 
aerosols (Dubovik et al., 2002) resulting in a higher absorp- 
tion coefficient. The main hypothesis then is that amongst 
the atmospheric aerosols measured at Neumayer dust is the 
only absorbing one. However, it should be noted that smoke, 
another absorbing aerosol, from South America has been ob- 
served in Antarctica (Pereira et al., 2006; Fiebig et al., 2009). 
This fact is considered in the analysis section below. 

2.2.3 Surface meteorological observations 

Weather meteorological stations located at several airports 
in Patagonia routinely monitor visibility conditions and stan- 
dard environmental variables in support of aircraft activity. 
These observations are carried out by operators who report 
a quantitative measure of sky turbidity (visibility , measured 
in kilometers or miles) and a code that identifies the state of 
the sky (such as the presence of precipitation, haze, smoke, 
etc.) . Dust events are explicitly reported. The data used 
to constrain the model simulations in this analysis is from 
the stations Comodoro Rivadavia (abbreviated CR and lo- 
cated at 45.8 S, 67.5 W, 46ma.s.l.) and Rio Grande (abbre- 
viated TdF and located in the island of Tierra del Fuego, 
53.8 S, 67.8 W, 22 m a.s.l). In addition, visibility data from 
the following stations was checked for dust activity: Bahia 
Blanca, Trelew, Puerto Deseado, San Julian and Rio Galle- 
gos. Visibility reports from meteorological stations consti- 
tute one of the longest record of modern dust activity and it 
can be used for dust trend studies (Mahowald et al., 2007). 
In this case, surface visibility information is used to obtain 
an independent confirmation of dust activity and the duration 
of the event. It should be noted that even though the sta- 
tion data is available NOAA’s National Climatic Data Cen- 
ter (http://www.ncdc.noaa.gov/oa/ncdc.html), the data used 
in this analysis was obtained directly from the argentine 


weather bureau (Servicio Meteorologico Nacional or SMN). 
The SMN data has the highest time resolution of observations 
(usually every hour) whereas not all the stations had the same 
time resolution at the NCDC archive. The use of higher tem- 
poral resolution made possible to resolve dust events of short 
duration. 

2.3 Hysplit model 

The HYSPLIT model (version 4.8) is a lagrangian trans- 
port model frequently used for trajectory analysis and dust 
and smoke forecasts (Draxler et al., 2001; Draxler and Rolf, 
2003). The reader is referred to Draxler and Hess (1998) for 
more details on the model and only those features relevant 
to this study are highlighted here. HYSPLIT uses the mete- 
orological fields from the NCEP’s Global Data Assimilation 
Scheme (GDAS) (Kanimatsu, 1989; Kalnay et al., 1996) at 
lxl degree resolution. Dust emission is based on the thresh- 
old friction velocity approach where surface roughness is 
correlated with soil properties. When the local wind speed 
exceeds the threshold velocity for the soil type of that emis- 
sion cell (Escudero et al., 2006), the dust emission rate is 
computed. The model domain was set to the southern hemi- 
sphere (0-90 S, 180 W to 180 E) at lxl degree resolution. In 
addition to the default dust emission parameters, the model 
allows the user to set a number of parameters such as size and 
intensity of the source and duration of emission. Model runs 
divided the atmosphere in 8 layers from 0 to 10 km. As ex- 
plained later in Sect. 3.1, all these options were set based on 
satellite and surface meteorological information. A few test 
runs determined that for a number of particles higher than 
2x 10 5 , the distribution of concentrations did not change and 
all runs used this number of released particles. A time step 
of computation of 12 min was used and output concentrations 
at surface and total in the column were saved at 1 hour inter- 
vals. The results are reported as number of particles per unit 
of volume in the layer or integrated over the column. 

It should be noted that model runs are used to indicate the 
location and distribution of the dust cloud as it propagates 
and to estimate the travel times to Antarctica. Because the 
unknown actual deposition rates, quantitative comparisons 
between model and in- situ/satellite retrievals will be avoided. 
However, some relative quantitative comparisons between 
the different runs will be shown to point out the strength of 
the sources and their relative contribution to each site. 

3 The modeling of dust events 
3.1 Identification of dust sources 

In order to determine the source regions, back- trajectories 
modeling runs were carried out from 9 March to 25 Febru- 
ary for air masses arriving to Dome C. Air masses were over 
the east coast of Patagonia, over the south tip of South Africa 
and over the south coast of Australia by 26 February. In order 
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to verify whether any of those areas were actively producing 
dust in the same time frame, the OMI AI, visible images and 
optical depths of MODIS were analyzed on these three areas. 
Three dust events were identified. One of them was in Cen- 
tral Australia near lake Eire (28.86 S, 139.13 E) on 2 March 
(see Aqua granule 03:05 UTC). This event was simulated by 
HYSPLIT in forward mode. The model run showed that air 
masses leaving Central Australia on 2 March advected east 
and did not reach Concordia. Thus, most likely Australia did 
not contribute to the dust measured at Concordia on 9 March. 

The other two dust events were observed in the Patagonia 
desert. Dust plumes were detected by MODIS on 26 Febru- 
ary 2005 (Fig. 2, Aqua granule 18:00) leaving the area of 
dry lakes between the city of Rio Grande (53.86 S, 67.76 W) 
and the San Sebastian Gulf in the island of Tierra del Fuego 
(about 90 km North of Rio Grande). The sources at the San 
Sebastian Gulf are deflation flats and are common in the area 
(Isla et al., 1991 ; Vilas et al., 1999) and the dry lakes north of 
Rio Grande are of the same nature (D. Gaiero, personal com- 
munication, 2009). A large plume is observed originating in 
a coastal mud flat or dry lake west of the San Sebastian gulf. 
The longest plume visible in the image originates from this 
source; its length is ~500 km and the optical depths are in the 
0.35-0.45 range, which are considered above “background” 
values (typically below than 0.15). Although not clear in this 
image, smaller dust plumes can be indentified north of the 
city of Rio Grande in the 250 m resolution image. This event 
lasted approximately 12 h based on the visibility reports from 
the meteorological station at Rio Grande. 

The second event occurred on 22 February and the dust 
originated from two different regions in Patagonia. Dust 
plumes were visible in MODIS (Aqua granule 18:30) im- 
agery originating at the shores of lake Colihue Huapi, 
(45.65 S, 68.75 W), west of the Comodoro Rivadavia sea- 
port. The source area has been identified as extremely windy 
(Fabraga, 1991) and as going through a severe desertification 
process (Del Valle et al., 1998). A compact patch of AODs 
(~0.3) was retrieved by MODIS over the ocean in the San 
Jorge Gulf (less than 100 km downwind from the source). 
The Tierra del Fuego island was cloud covered in the Aqua 
and Terra-MODIS images but the sea along the shore was 
cloudless. The MODIS aerosol product for this clear sky 
area showed a large group of pixels of high and homogenous 
aerosol optical (AOD~>0.3, Aqua granule 18:25) approxi- 
mately 100 km SE of the San Sebastian Gulf. Dust activity 
was confirmed by surface weather stations in both Comodoro 
Rivadavia and Rio Grande. 

No biomass burning or additional dust activity was noted 
in the 22 to 26 February period in the satellite imagery and 
in the surface meteorological observations. 

3.2 Setup of forward model runs 

Forward model runs with HYSPFIT were simulated for the 
22 and 26 February dust events. The runs were initialized 



Fig. 2. Dust Event in the Tierra del Fuego Island (~53 S) detected 
by MODIS in the Aqua satellite on 26 February 2005. The largest 
plume is in the Gulf of San Sebastian and several minor plumes 
(some are not resolved in this image but they are visible by zooming 
in the full 250 m image) are located between San Sebastian and the 
city of Rio Grande, 100 km south. 


with the information from the known (observed) dust sources 
and the duration of event. The exact location and the area 
of the dust sources were digitized from the MODIS images. 
High-resolution MODIS images (250 m) were used to locate 
each dust emitting lake. The source was identified using 
Google Earth (GE) and then a polygon was drawn over the 
identified sources. If the MODIS image was clear enough, 
the section of the lake emitting dust was drawn (for exam- 
ple, a shore of a lake). Otherwise, the whole area of the lake 
was drawn and it was assumed that the whole lake was ac- 
tive. Then, the area of the polygon was computed and used 
as one of the inputs for initialization of the model simula- 
tions. The advantage of this approach resides in model runs 
based on identified active sources and the length of time of 
emission is based on the weather station reports of observed 
dust. There is, of course, an error in estimating the active 
area in the GE imagery based on the dust activity seen in the 
MODIS image. For example, the edges of the sources in the 
GE image database may not correspond to the edges of the 
active sources in MODIS. These differences are assumed to 
be small. Table 1 details the location of each source identi- 
fied in the MODIS Terra and Aqua image and their respective 
area in square kilometers. 

The duration of each event was set by the number of hours 
where any amount of dust was observed (i.e. World Mete- 
orological Organization METOP codes 6-9) by the closest 
meteorological station to the source. None of the stations 
reported a “dust storm” (code 9) and most corresponded 
to codes 6 and 7 (low and moderate dust in the air). For 
each group of active sources identified, a HYSPFIT run 
in forward mode was carried out. These runs are identi- 
fied as Feb22TdF, Feb26TdF and Feb22CR corresponding 
to dust leaving Tierra del Fuego on 22 and 26 February and 
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Table 1 . Location and description of dust activity and sources as seen in MODIS visible images (observation time noted as “granules”). In 
all cases when it was not clear if dust is coming from the edge of the source, the estimated area of emission was to be the totality of the 
lake/depression. 


Date 

22/Feb/2005 

Terra granule 14:20 UT C 

Latitude 

Longitude 

Area (km 2 ) 

Remarks 

-45.66 

-68.66 

40.194 

Patagonia is located in center of MODIS swath, providing good detail. Area is generally cloudy, 

-53.23 

-68.68 

15.7180 

with a clear area extending from Comodoro Rivadavia to Puerto Santa Cruz. North end of 

-53.09 

-68.58 

2.2161 

Tierra del Fuego is partly cloudy. Small lakes appear to have minimal dust activity. 

-53.63 

-68.15 

6.8615 


Date 

22/Feb/2005 

Aqua granule 18:25 and 18:30 UTC 

Latitude 

Longitude 

Area (km 2 ) 

Remarks 

-45.66 

-68.66 

40.194 

18:30 UTC: Clear dust cloud off the coast of Comodoro Rivadavia. Some dust activity from 
same source seen in Terra. 18:25 UTC: Area of interest is located on the far western edge 
of image and it is mostly obscured by patchy clouds south of San Jorge Gulf. Clear sky off 
the coast of Tierra del Fuego. Dust off the coast of Rfo Grande. 

Date 

26/Feb/2005 

Terra granule 1 4 : 00 UT C 

Latitude 

Longitude 

Area (km 2 ) 

Remarks 

-53.23 

-68.68 

15.7180 

Clear sky view over Tierra del Fuego but towards the west edge of then MODIS swath, causing 

-53.66 

-68.26 

7.5952 

slight loss of detail. 

-53.78 

-67.79 

2.3015 


-53.86 

-67.62 

9.2771 


-53.18 

-69.11 

9.4524 


-53.58 

-68.34 

4.2668 


Date 

26/Feb/2005 
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Comodoro Rivadavia on 26 February. The duration of emis- 
sion is set to the time frame of dust activity as reported by the 
surface observations. These runs are referred as “constrained 
runs”. In order to test the impact of initializing the runs 
with observations, three additional runs for the same group 
of sources was carried out, with the default dust emission 
enabled i.e. sources become active when the threshold wind 
speed is reached any time between at 00:00 and 24:00 UTC 
on 22 and 26 February. These runs are referred as “uncon- 
strained runs”. Thus a total of six model runs were carried 
out. 

According to the weather observations in Rio Grande 
airport on 26 February, dust was observed from 13 (local 
time = UTC-4 h) until 24 h (local time) with visibilities rang- 
ing 2-8 km in the first 8 h. The weather code assigned to this 
event was 7 which is “dust in the air”, that is, it was not iden- 
tified as dust storm. Wind direction was from the west and 


averaged 40 to 50 km/h and gusts reached 75 km/h. Then, the 
constrained simulations of HYSPLIT used a window of 11 h 
for active dust emission. Note that it is unknown for how 
long the largest source at San Sebastian was active since the 
weather station was too far to see this plume and visibility 
reports were most likely based on dust activity in local lakes. 
However, it was active between 14:00 and 19:40 UTC since 
dust emission was clear in the MODIS images. 

As opposed to 26 February, most of Tierra del Fuego was 
cloud covered on 22 February during both Terra and Aqua 
passes. However, the MODIS aerosol retrieval algorithm de- 
tected an aerosol cloud in the cloudless area near the shore of 
Tierra del Fuego. In addition, the Rio Grande meteorologi- 
cal station reported reduced visibility by dust from 14 to 19 h 
local time with wind speeds ranging from 50-60 km/h and 
wind gusts ranging 75 to 95 km/h. Since it was cloudy, there 
is no visual imagery to confirm the location of active sources. 
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Fig. 3. Total columnar dust concentrations modeled by HYSPLIT 
initialized by the active dust sources detected by MODIS on 26 
February 2005 and emission duration constrained by surface visi- 
bility measurements. Outputs are shown for 27 February (left), 3 
March (center) and 9 March (right). Lighter colors indicate higher 
concentrations. 

Thus the model runs assumed the same active sources visible 
in a previous dust event in the same area (16 January) which 
are very close in size and proximity to those seen active on 26 
February. The Colihue Huapi lake area was clear of clouds 
in the MODIS Aqua image (granule 18:30UTC) and there 
was positive identification of dust emission from the eastern 
shores of the lake. The Comodoro Rivadavia weather station, 
located just downwind of the Colihue Huapi lake, reported 
dust activity from 13 to 20 h local with gusts in the 65 to 
85 km/h range during the period. Overall and according to 
surface observations, it appears the dust event of 26 Febru- 
ary was the most intense and the longest amongst the three 
events. 

Although in the reminder of the paper the constrained runs 
will be used for comparison with satellite and surface data, 
there were some differences between the two modes of emis- 
sion that will be noted. 

3.3 Propagation of the modeled dust cloud 

Figure 3 shows normalized total columnar concentrations for 
the case Feb26TdF for 2, 6 and 11 days after the event. The 
rapid progression is apparent on the second day after emis- 
sion with the leading edge located +1000 km from the source. 
By the 5th day, the dust cloud had reached over the area be- 
tween South Africa and Antarctica and the leading east edge 
had reached Kerguelen Islands (70 E) area. The south end of 
the cloud is entering Antarctica at ~20E. By the 11th day, 
the dust covered most of the Antarctica plateau. It should be 
noted that after the 5th day, the dust has dispersed and di- 
luted significantly and the area covered is very large. This 
is not necessarily is a proof of actual dust transport to that 
point since the removal processes are poorly parameterized 
in models (see for example Xian et al., 2009) and it is un- 
known how much material removed from the column. How- 
ever, the model output does provide an envelope of maximum 
expansion of the dust cloud and a minimum travel times un- 
der the ideal condition of no removal events. 
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Fig. 4. Time of HYSPLIT surface concentrations at Concordia 
(123 E, 75 S) for 3 dust events described in text and the dashed black 
line is the nss-Ca filter measurement at Dome C. Top: the contribu- 
tion of each run to the total concentration for that day and Bottom: 
Surface concentrations normalized to the maximum of Feb26TdF. 
Top figure illustrates the relative contribution of each dust event. 
Bottom figure shows that most of the dust arriving to Dome C on 8 
March originated in Tierra del Fuego on 26 February. 


3.4 Comparison with surface measurements 

A comparison of model outputs with surface measurements 
are shown in Fig. 4a-b where a time series of the dust con- 
centrations from the lowest layer at the Concordia site for the 
three constrained simulations is displayed. Figure 4a shows 
the contribution of each source to the total dust concentration 
at the site (modeled data is aggregated to 24 h intervals). Fig- 
ure 4b shows all runs normalized to the maximum value of 
the Feb26Td run. The dashed line is the normalized surface 
measurement of nss-Ca at Concordia. 

Figure 4a shows that dust from event Feb22TdF starts to 
arrive on 3 March and for event Feb22CR arrives in full on 6 
March. The different travel time of the two sources (~9 days 
and ~12 days respectively) is a consequence of the different 
location of the sources (+800 km apart). In particular, dust 
from Comodoro Rivadavia has to cross the storm track where 
emission in Tierra del Fuego is already emitted at the track’s 
path or south of it. It takes ~9 days of travel time for the 
event Feb26TdF too. 

Most notably, the maximum amount of dust arriving for 
run Feb26TdF occurs during most of the day on 8 March 
that it is, during the 24 h filter exposure at the site. The 
Feb26TdF event contributed much more dust to the Concor- 
dia site than any of the other two events. Figure 4a shows 
that the 22 February events contributed most of the dust ar- 
riving to Concordia prior to 8 March but the amount con- 
tributed by this event appear to be minimal when compared 
to the other two events. According to the model, the concen- 
trations from Feb22TdF were always 2-3 or more orders of 
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Fig. 5. Time series of HYSPLIT surface concentrations at 
Neumeyer for the 3 simulated dust events described in text and the 
black line is the normalized aerosol absorption coefficient measure- 
ment at the site. Top plot: the contribution of each run to the total 
concentration for that day and bottom: concentrations normalized to 
the maximum of Feb26TdF. The 26 February event correlates with 
increases in absorption on 5 and 9 March. Dust from the Feb22CR 
event arrives earlier to Nuemayer (1 March). 


magnitude than the other two cases (this is the reason why 
there is no visible blue line in Fig. 4b). 

Figure 5 shows a similar time series for the Neumayer sta- 
tion. In this case, comparisons are not as clear cut. Aerosol 
absorption observations (black dashed line) are normalized 
to the 9 March value and indicate the presence of three peaks 
(28 February-1 March, 4-5 and 9-10 March). According 
to the model runs, dust from two sources (Feb22CR and 
Feb26TdF) arrive on 4 and 9 March and they coincide with 
an increase in aerosol absorption (Fig. 5a). However, note 
that model runs show dust from Feb22CR and Feb26TdF ar- 
riving on 7 March too but no corresponding increase in ab- 
sorption is observed and there is no coincidence with the 28 
February-1 March peaks either. Figure 5b shows the same 
model data but normalized to the maximum value of the 
Feb26TdF for the corresponding day. After 3 March, dust 
from Feb26TdF is the dominant contributor to Neumayer 
and CR was a minor contributor. Like in the previous com- 
parison, the Feb22TdF contribution with respect to the other 
larger events is minimal. In addition, the lack of coincide be- 
tween some of the peaks suggests that additional absorbing 
aerosols other than dust reached Neumayer in this period. 

This travel time may appear long if considering that the 
Neumayer is relatively close to Patagonia (compared to Con- 
cordia). However, inspection of the modeled trajectories 
for each of the three cases shows that none of them fol- 
lows a straight meridional transport. This transport pattern is 
common in the high latitude SH where moisture fluxes into 
Antarctica follow the motion of cyclonic systems (Jones and 
Simmonds, 1993; Rassmussen and Turner, 2003). In all three 
cases the dust cloud is advected east and then it veered south 


and then west arriving to Neumayer from the NE (see for 
example Fig. 3, center). This path results in a longer trans- 
port time than would be expected if going in a more or less 
straight line from Patagonia to Neumayer. This transport pat- 
tern is partially confirmed by the satellite imagery (see next 
section). 

These two comparisons suggest that dust leaving Tierra 
del Fuego on 26 February contributed most of the dust ob- 
served in Dome C on 8 March and it was the largest contrib- 
utor to the aerosol absorption measured in Neumayer after 3 
March. Although the sources at Comodoro Rivadavia have 
larger emission area than Tierra del Fuego, it appears the 
proximity to Antarctica makes the latter region a more ef- 
fective source of dust in quantitative terms. The travel times 
to Neumayer for the Feb22TdF and Feb22CR and Feb26TdF 
runs are ~7, ~11 and days respectively. 

Figure 4a illustrates another interesting point. If signifi- 
cant amounts of dust leaving Tierra del Fuego on 26 February 
reached Antarctica, why did not the case of Feb22TdF show a 
similar pattern? Since both runs had sources of similar size, 
the other two controlling parameters of dust emitted in the 
model are wind speed and the length of time exceeding the 
threshold velocity. The amount of dust generated is propor- 
tional to a power function of the wind and small differences 
in wind speeds have a non-linear effect and large variations in 
dust production. Assimilated wind fields and measured sur- 
face winds at the Rio Grande station agree in showing more 
intense winds and longer emission times on 26 February than 
on the 22nd that is, the model produced more quantities of 
dust on 26 February. Thus, under the assumption that in both 
events the removal rates of dust were similar during the tran- 
sit to Concordia, the number of hours of wind speeds above 
the threshold velocity are the critical parameters in making 
the Feb26TdF event the major dust contributor. 

The comparison of assimilated wind speeds (GDAS data 
base) in HYSPLIT for the grid point corresponding to Rio 
Grande station with the measured winds at the same site 
shows that high winds start earlier in the assimilated winds 
than at the surface station. This time difference in the initi- 
ation of high winds impacts the dust production of the un- 
constrained and constrained model runs. The unconstrained 
runs emit larger quantities of dust and earlier than in the con- 
strained runs. These differences do not appear to have an im- 
pact on the travel times to Neumayer and Concordia. How- 
ever, the differences become obvious when comparing with 
satellite as it is shown in the next section. 


4 The satellite record 

The previous section showed that dust leaving Patagonia can 
reach sites in Antarctica more than 5000 km from the source. 
However, the long time difference between emission and ar- 
rival to the Antarctica sites as well as the poor characteriza- 
tion of the removal processes commonly found in transport 
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Fig. 6. Hemispherical map aerosol optical depths from Terra MODIS on 26 February 2005. Note the adjustment of color scale by setting 
maximun at 0.3 in order to emphasize the variable low aerosol concentrations. The dust event from Tierra del Fuego is visible with AOTs 
>0.24 (see text). 


models casts some doubts on whether the dust from Patago- 
nia c an reach this area. In this section the analysis of satellite 
observations of the dust transport from Tierra del Fuego (26 
February event) is described in detail. 

4.1 Methods 

The abundant cloudiness over the Southern Ocean (SO) is 
a significant shortcoming because aerosol remote sensing re- 
quires a clear sky patch (i.e. a cloudless view of the ocean) for 
the retrievals algorithm to operate. This is illustrated in the 
MODIS-Terra AOTs (Fig. 6) where large areas in the SO are 
devoid of aerosol retrievals because of clouds. Also, note that 
in order to illustrate the presence of aerosols the color scale 
has to be set from 0.05 to 0.3, a range of values normally con- 
sidered clear to moderate concentrations of aerosols. Further, 
this image includes the initial stages of the 26 February dust 
event in Tierra del Fuego. However, its AOT values are in the 
0.2-0. 3 range and they are hardly any different to the equally 
high values of AOT found in the rest of the scene. With- 
out any additional information, it is not possible to identify 
a plume or cloud of dust in this image. The most common 
aerosol types present in this region are commonly referred as 
“marine” aerosols. The optical properties of marine aerosols 
are typical of those aerosols dominated by a coarse mode 
resulting in low spectral dependence in the visible and low 
aerosol optical depths (typically below 0.15). Dust detec- 
tion approaches (Kaufman et al., 2005; Yu et al., 2009) dif- 
ferentiate dust from seasalt by using a threshold based on 
the AOT magnitude. This is a reasonable assumption in en- 
vironments where the abundant and constant dust emission 
results in high AOTs (>0.4). However, they are not well 
suited to study low aerosol concentrations and heavy cloudy 
conditions such as those encountered in the transit to Antarc- 
tica. Thus, Fig. 6 illustrates one of the existing difficulties 
of aerosol detection and identification in the Southern Ocean 
using visible-infrared detectors such as MODIS. In this case, 


the high spatial resolution of the pixels used by the MODIS 
algorithm (250 m channel, Martins et al., 2002) performs rea- 
sonably in detecting clear sky patches. Yet, the aerosol con- 
centrations are not high enough to warrant a method such as 
the one used by Kaufman et al. (2005) to identify aerosol 
types. Additional information regarding aerosol identifica- 
tion is needed. 

A number of satellite detectors were available at the time 
of the dust events here studied. In addition to MODIS, we 
explored the level 2 aerosol retrievals of the following de- 
tectors: OMI, PARASOL, MISR and GLAS. Although OMI 
and PARASOL have been used for aerosol type identifica- 
tion, their retrievals during the period 22 February to 9 March 
in the Southern Ocean region did not show useful informa- 
tion for this analysis. The MISR detector on-board of Terra 
has a narrower swath than OMI and PARASOL. Level 2 
Aerosol retrievals are reported at M 7 x l 7 km but the re- 
trieval is carried by aggregating smaller pixels (170x170 m) 
and selecting only those where the clear sky confidence is 
high (Kahn et al., 2005). Thus, Level 2 Aerosol retrievals 
of non-spherical aerosol fraction are representative of pix- 
els much smaller than those used by OMI and PARASOL. 
In addition, MISR data is available on the same days and 
within the same swath of MODIS-Terra. All these features 
made MISR the tool of choice for aerosol identification in 
this study. 

Although operating a couple of months per year, the 
GLAS lidar was active during late February and March 2005 
and it provided the vertical distribution of aerosols. Because 
operational problems, the 532 nm laser was not working 
properly and did not provide useful profiles. The 1064 nm 
laser was operational and it was sensitive to aerosols in clear 
sky conditions. However, the laser was losing power too 
and no aerosol layers were identified above clouds in part 
because the contamination of the return signal from clouds 
(Steve Pal., GLAS science team member, personal commu- 
nication, 2009). 
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Fig. 7. MODIS-Terra AOTs (left panel) and MISR non-spherical fraction (right) for pixels with AOTs >0.16 (pass time near 13:55 UTC). 
The MODIS image includes the envelope of pixels witn non-zero dust concentration according to HYSPLIT for a model run unconstrained 
by surface observations (thick dashed line). Also, the image shows the north- south track of GLAS five hours later. The image on the right 
shows the envelope of non-zero dust concentrations (solid thick line) for HYPSLIT model run constrained by surface obs. 



1 .0e-05 
3.2e-06 
I 0e-06 
3.2e-07 
1 0e-07 


HYSPLIT output averaged over 3 hours. Run Feb26Tdf 



Fig. 8. Top Panel: GLAS 1064nm backscatter profiles for the low 6.6 km. White=no signal, red=clouds, yellow and green; orange=aerosol. 
Bottom panel: normalized HYSPLIT concentrations for run TdF26, model run constrained by observations. 


The next sections show MODIS, MISR and GLAS (when 
available) observations analyzed jointly with the model out- 
puts for run Feb26TdF. 

4.2 Analysis of case 26 February 

Figures 7 show retrievals for 26 February. Aerosol optical 
depths from MODIS-Terra are in the left figure and MISR 
non-spherical AOD fraction on the right. Because MISR is 
onboard of the Terra satellite, both MODIS and MISR obser- 
vations are simultaneous. MODIS data is provided in 5 min- 
utes segments (“granules”) and this image shows the aggre- 
gation of two consecutive overpasses or orbits. Data of each 
overpass corresponds to observations at the same local time. 
The MISR data is only displayed for those pixels with AOTs 
above 0.16. Non-spherical retrievals below this threshold are 
considered unreliable (Kalishnikova and Kahn, 2008). Over- 
lapped to both images is the contour of the dust cloud mod- 
eled by HYSPLIT (averaged over the 3 hours closest to the 


satellite pass). The right hand side shows the contour for the 
constrained run and the unconstrained run is on the left. The 
contour displays all pixels with non-zero dust concentration 
according to the model. The image in the left includes the 
track of GLAS which flew over at approximately 19:50UTC. 
The GLAS vertical profiles of backscattering and the re- 
spective HSYPLIT output profiles along the GLAS track are 
shown in Fig. 8. Qualitatively, aerosols can be distinguished 
by the colors green, yellow and orange colored (~8.10 -7 to 
3.10 -6 1/m.str) and clouds are in red. White colors are neg- 
ative values or no signal return from under the clouds. 

By comparing both MODIS and MISR retrievals, there is 
a clear spatial consistency between the high AOTs and the 
high values of non-spherical fraction. In addition the loca- 
tion of the dust cloud according to HYSPLIT coincides with 
the spatial distribution of the aerosol according to both satel- 
lites. However, there are notable differences between the 
constrained and unconstrained model runs. The contour of 
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Fig. 9. Same as Fig. 7 but for the day after of the dust event (27 February). GLAS pass (N-S dashed line) was at ~18:20UTC. The MODIS 
images are composed out of two successive orbits in the region (~ 11:20 and 13:00 UTC). 
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Fig. 10. Same as Fig. 8 but for 27 February. Note that the location 
does not coincide with location of the dust modeled by HYSPLIT. 

the unconstrained run seems to contain most of the MODIS 
higher optical depths whereas it is not the case for the con- 
strained run. By inspecting and comparing the assimilated 
data used in the unconstrained run (GDAS) with the surface 
observations, the latter report visible dust and high winds for 
fewer hours than the number of hours of high winds in the 
assimilated data. In other words, the unconstrained run starts 
to emit dust earlier and for a longer length of time than the 
constrained run. 

The GLAS profile shows high backscattering (orange and 
yellow) in the bottom 2 km of the profile and it extends from 
49 to 53S consistent with MODIS and MISR spatial distribu- 
tion of the aerosol in the north- south direction. The vertical 
location of the dust is consistent with the model simulation 
that places the dust in the bottom two layers of the model 
(bottom Fig. 8). The MODIS-Aqua pass was closer in time 
to the GLAS overpass and had similar AOTs and aerosol dis- 
tribution as derived by MODIS -Terra (not shown). 

In summary, on the day of the event MISR, MODIS, 
GLAS and the model simulation are all consistent in plac- 
ing the dust in same region and height. This agreement is 
useful for the analysis of satellite data in the subsequent days 
where dust concentrations are lower and the dust cloud is 


of highest backscattering seen by GLAS in the boundary layer (~44 S) 


more diffuse and difficult to identify just by looking at the 
spatial pattern. 

4.3 Analysis of the day after the event (27 February) 

Figures 9 and 10 show a set of MODIS, MISR and GLAS 
images in similar format to Figs. 7 and 8 for the day after 
of the event (27 February). In this case, the dust cloud has 
dispersed and it cannot be easily identified in a visible im- 
age (not shown) like the day before. Thus all observations 
need to be considered jointly in order to detect and locate 
the extension of the dust cloud. In this case, Fig. 9a and b 
and model contour are consistent too. The MODIS AOTs 
are higher inside the model contours than the surroundings. 
In the same area, MISR shows several pixels with an impor- 
tant non- spherical fraction contribution in the total optical 
depth (ranging 38-74%). In addition, MISR retrievals are in 
agreement with the location of the dust cloud according to 
the model. 

The GLAS image (Fig. 10a) shows high aerosol concen- 
trations in the boundary layer as far south as 53 S. The model 
output places the dust in the boundary layer too in the same 
area (see the track in Fig. 9a). Further, the north to south 
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Fig. 11. MODIS-Terra AOTs (left) and MISR non- spherical fraction (right) for 3 March. MODIS images are composites of two consecutive 
overpasses (at 10:55 and 12:35 UTC). North to South dashed line is the track of GLASS at 14:13 UTC. 


extent of the dust cloud agrees with GLAS too although the 
model does not extend as far north as the GLAS signal sug- 
gests. There are some inconsistencies too, note the high 
backscattering at 44.5 S does not seem to be dust according 
to the model. 

As noted in the previous sections, the difference in emis- 
sion schemes in the HYSPLIT runs results in different spa- 
tial distribution and reach of the dust cloud. In this case, 
both constrained and unconstrained runs agree in the loca- 
tion and distribution of dust in the region enclosed by the 
non-spherical retrievals of MISR (50 W to 45 W). However, 
both runs disagree on the location of the leading (or east) 
edge of the dust cloud. The difference of emission initiation 
is more apparent after 24 hours. The MODIS image show a 
large patch of AOTs ranging 0.16-0.22 in the 38-32 W, 45- 
49 S box which agrees with the unconstrained run but it does 
not agree with the constrained run. This is a demonstration of 
the effect of the difference in emission conditions as defined 
by the constrained and unconstrained runs. 

When putting all this information together i.e. agreement 
between model height and location of the dust with aerosol 
information from the satellites, it becomes apparent that the 
dust cloud reached the south-central South Atlantic extend- 
ing from ~45 S to 50 S and reaching as far east as ~40W or 
further east. This agrees with an observation of dust previ- 
ously reported in the same area (Gasso and Stein, 2007). 

4.4 Analysis of days after 27 February 

By 28 February, the simulated dust cloud has extended 
throughout most of the South Atlantic. However, because 
the extensive cloudiness, there were very few MODIS and 
MISR retrievals within. Even when clear sky patches were 
found by the respective algorithms, they contained small or 
zero non-spherical fractions and low AOTS (0.05-0.1 range). 
On 1 and 2 March, there were large clear sky patches within 
the simulated dust cloud with high AOTs (0.15-0.25). How- 
ever, none of them contained a high number values of non- 
spherical fraction. After 4 March, analysis of images and 
retrievals of subsequent days did not show any positive iden- 
tification of dust. 


On 3 March, however, two large patches of clear sky 
within the (constrained simulation) dust cloud were located 
SW of South Africa (35-55 S, 25W-20E). Optical depths 
within both patches are well above 0.1 with several pix- 
els above 0.3 (Fig. 11). Further, each of the patches has 
dozens of contiguous clear sky pixels with a gradual change 
of AOT, which is a good visual indication of minimal cloud 
contamination. There were also abundant retrievals of high 
values of non-spherical fraction (ranging 0. 5-0.7) in those 
clear sky patches. Further, four hours later Aqua (near 
13:25 UTC) observed similarly high AOTs in the same two 
patches (Fig. 12a). Approximately at 14:13 UTC GLAS flew 
by over the same clear patch observed by Aqua (see track 
overlap in Fig. 11). The GLAS profile confirms of the pres- 
ence of abundant aerosol concentration in the boundary layer. 

In order to get additional insight of the nature of the high 
concentration of non-spherical aerosols, a few additional 
checks were performed. Because cirrus are a common occur- 
rence and can interfere in the MISR aerosol retrieval, their 
presence in the area was tested. An image of the MODIS 
(band 26) reflectance at 1.64 pm MODIS was generated. This 
band is particularly sensitive to high altitude cirrus (Acker- 
man et al., 1998) and to water vapor absorption. If cirrus are 
present, they scatter back radiation to the sensor. However, 
if there are no cirrus, the signal is absorbed by the boundary 
layer water vapor resulting in an almost zero reflectance for 
this pixel. The corresponding images clearly showed very 
low and almost zero band 26 reflectances in the clear sky 
patches where the high AOT were retrieved. This suggests no 
cirrus contamination. The large size of the patches observed 
on this day prompted the reexamination of the OMI images 
which yielded a negative result (i.e. no AI>1 was detected 
in the area of high values of non-spherical fraction). This is 
not entirely inconsistent because OMI is not very sensitive to 
low altitude dust (Mahowald and Dufresne, 2004), which is 
the case here as the GLAS data for this day indicate. 

An additional test was applied using a dust detection 
method based on IR emissivities (Ackerman, 1997). Because 
dust has a different spectral emissivities in the 11 and 12 pm 
(Bands 31 and 32 of MODIS), several detection techniques 
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Fig. 12. Top: MODISAqua Aerosol optical depths for 3 March (around 13:25 UTC) with the contour of dust concentration of constrained 
run Feb26TdF. The dashed line is the GLAS track for an overpass pass at 14:13 UTC. Bottom: the backscattering profile for the track shown 
in the AOT image and shows that the high optical depths correspond to boundary level aerosols 


have been developed based on these IR channels in combi- 
nation with visible channels (Baddock et al., 2009). In this 
case only the Ackerman test was applied. The test yielded 
no identification of dust in the infrared. However, the nega- 
tive result of this test does not rule out the presence of dust. 
Since in this case, the aerosol is located in the low 700m of 
the boundary layer (confirmed by GLAS) and the contrast in 
brightness temperatures may not be enough. Further there 
have been reports of brightness temperature differences op- 
posite to what the Ackerman et al., 1997 method predicts 
(Darmenov and Sokolik, 2005). 

In summary, the model simulation and MISR non- 
spherical fraction are consistent in pointing to the presence 
of dust in this area. MODIS confirms the presence of high 
concentration of aerosols and they are mostly restricted to 
the boundary layer. However, other dust detection proxies do 
not positively associate the aerosol with dust although these 
tests were not designed for this type of scene (dust over ocean 
in the boundary layer). 


5 Discussion 

Although ideally more sophisticated techniques (such as iso- 
topic analysis) would be able to trace the sources of dust 
measured at Concordia, there are a number of observations 
and modeling features that show a remarkable consistency in 
showing dust transport to Antarctica. The findings described 
in the previous sections can be summarized as follows: 

1 . Identified dust in the Concordia site on 8-9 March 2005 
(Fig. 1) 

2. Backtrajectory modeling indicates Patagonia as a possi- 
ble source (Fig. 3). 

3. Surface visibility and satellite confirm dust activity in 
Patagonia (Sect. 3.1) coinciding with the departure of 
air parcels that arrived to Concordia (backtrajectory 
analysis) 

4. Surface measurements of aerosol absorption in Neu- 
mayer show increases on 28 February to 5 and 9 March 
(Fig. 5). 

5. Forward modeling initialized by observed dust sources 
shows agreement in the timing of arrival of modeled 
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dust with the observed increase of dust observations and 
proxies (model comparison with surface data at Neu- 
mayer and Concordia, Sect. 3.4). 

6. Satellite observations confirms the presence of 
dust (traveling in the boundary layer) at least to 
~ 1700 km east from the source and possibly reaching 
~1800 km SW of Cape Town, South Africa (Sect. 4). 

The interpretation of the different agreements between 
model and observation does have some weaknesses. Some of 
tools used here have not been validated. For example, the un- 
certainties in assimilated winds and in deposition of aerosols 
have not been properly characterized in HYSPLIT. Thus, be- 
cause the lack of information of removal processes during 
transport, the approach used was limited to temporal com- 
parisons (i.e. transit times to Neumayer and Concordia) and 
quantitative comparisons were avoided when comparing with 
surface data. Only relative comparisons between normalized 
concentrations of different model runs were examined. How- 
ever, in comparing modeled concentration with satellite data, 
the criterion was to display all areas where dust could have 
arrived according to the model. Thus the model provided the 
envelope of maximum extension of the dust cloud and there 
was no attempt to compare the modeled concentration mag- 
nitudes with the satellite data. This assumption was reason- 
able for the first two days of the event when unknown errors 
in emission details and winds have not propagated. However, 
it is not clear the impact of these errors after 5-6 days after 
emission. An example of the uncertainty in emission con- 
ditions is demonstrated in the differences in the propagation 
and spatial distribution of the dust clouds in constrained and 
unconstrained model runs (Figs. 7 and 9). 

Also, when comparing model outputs with satellite data, 
the unconstrained model runs compared better with satellite 
data than the constrained runs. This suggests that surface 
information can have limited use unless the visibility obser- 
vations are close to the major sources, unlike the case studied 
here. 

Satellite data has its limitations too. Although aerosol op- 
tical depth (particularly from MODIS) is a parameter exten- 
sively validated in regimes of biomass burning, dust and pol- 
lution aerosols, they have not been validated as extensively 
in the marine environment in particular at high latitudes. In 
fact, in the case of MODIS, some biases of AOT retrievals in 
the marine environment have been noted (Zhang and Reid, 
2005; Zhang et al., 2006). Because MODIS data was used 
in conjunction with MISR retrievals (which uses a different 
approach than the MODIS retrievals) and in all cases shown 
here, they showed agreement indicating no particular bias in 
MODIS AOT. Finally, the MISR non-spherical fraction pa- 
rameter has not been validated. However, the MISR team has 
used this parameter to characterize long range transport of 
Saharan dust (Kalashnikova and Kahn, 2008) with the same 
approach that was used here. 


As it is apparent in Figs. 6,7 and 9, there are several clus- 
ters of AOTs in the same range of values as those identified 
as dust. Likewise with the MISR non-spherical fraction im- 
age, there are several spots of high values outside the model 
contours similar to those associated with the dust cloud. The 
identification of the dust becomes only apparent when all 
these pieces of information (model and satellite) are analyzed 
together. Specifically, it is critical in this analysis the identifi- 
cation of a dust event at its initiation, which allows the correct 
setup of the model run and obtains an accurate prediction of 
the location of the dust cloud. 

When using OMI data, the rule of thumb for dust detection 
is AI>1 and values in the range 0-1 are normally dismissed 
as not reliable as a dust flag. Reexamination of the OMI Al 
retrievals for the time period 27 February to 4 March shows 
an interesting feature. In all clear sky patches where high 
non-spherical fraction was detected, the corresponding Al 
for approximately the same area (given the time difference 
between MISR and OMI passes) was in the 0 to 1 range. 
This suggests that there was an absorbing agent either in the 
ocean or the atmosphere. If it was in the atmosphere, it was 
a low laying and/or low concentration of absorbing aerosols. 
If it was originating in the ocean, it could be water absorp- 
tion or absorption by biogenic materials (Omar Torres, OMI 
Aerosol team member, personal communication, 2009). Al- 
though an interesting possibility for further study, this is an 
analysis that is beyond this work. 

The GLAS data provided useful information on the loca- 
tion of the aerosol load in the atmospheric column. A sig- 
nificant amount of dust remains in the boundary layer and 
for a long distance and it confirms the modeling work of Li 
et al, (2008) and Krinner et al, (in press) who found signifi- 
cant transport of dust occurs in the boundary layer. It should 
be noted, however, because deteriorating laser conditions in 
GLAS, significant number of profiles in cloudy conditions 
could not be used. As a result, it is not possible to rule out 
dust transport above clouds in the free troposphere. Although 
it is clear that the dust event studied here remained in the 
boundary layer in the first two days, there are no lidar ob- 
servations south of ~55S that can corroborate the location 
of the dust in the column. Our model simulation suggests the 
dust is lifted as it moves south but there is not observable that 
can verify this. With this regards it is clear that further case 
studies need to be studied in the CALIPSO (a spacebome li- 
dar deployed in 2006 i.e. after the dust events studied here) 
data set. 


6 Conclusions 

The purpose of the paper is to show the difficulty to char- 
acterize dust over the Southern Ocean. This effort com- 
prised the joint analysis of surface measurements of aerosols, 
aerosol dispersion modeling and characterization and track- 
ing of the dust cloud in its transit over the Southern Ocean. 
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For the case studied, it is clear that dust emission can occur 
as far south as latitude 54 S, making the Tierra del Fuego Is- 
land the southernmost dust source currently active. Also, it 
appears that dust from Patagonia reached Antarctica in early 
March 2005. This is supported by the number of pieces of 
evidence (listed in Sect. 5) that point to long range transport 
between the two continents. 

Even the combination of information from several satel- 
lite observations (many of them considered the best tools 
for aerosol studies currently) do not permit an unambigu- 
ous determination of the aerosol type under observation. As 
noted through the paper, there are numerous ambiguities in 
the satellite and modeled data that can yield mixed results in 
terms of characterizing the transport. Only when modeled 
and satellite data are interpreted jointly, a more consistent 
description of the dust cloud becomes apparent. 

Also, it is clear from this study that dust from the Patag- 
onia desert reach iron deficient regions in the SO and that 
much of this dust remains in the boundary layer suggesting 
the bulk of the dust cloud deposits over the ocean surface. 

Finally, this works shows that modern dust activity in 
Patagonia can reach Antarctica. Thus, by studying current 
dust transport in the SO, better insights can be gained on the 
understanding of dynamics and dust deposition in past glacial 
periods. 
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